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We study special relativistic effects on the entanglement between either spins or momenta of com- 
posite quantum systems of two spin-| massive particles, either indistinguishable or distinguishable, 
in inertial reference frames in relative motion. For the case of indistinguishable particles, we consider 
a balanced scenario where the momenta of the pair are well-defined but not maximally entangled 
in the rest frame while the spins of the pair are described by a one-parameter (rf) family of entan- 
gled bipartite states. For the case of distinguishable particles, we consider an unbalanced scenario 
where the momenta of the pair are well-defined and maximally entangled in the rest frame while 
the spins of the pair are described by a one-parameter (f ) family of non-maximally entangled bipar- 
tite states. In both cases, we show that neither the spin-spin (ss) nor the momentum-momentum 
(mm) entanglements quantified by means of Wootters' concurrence are Lorentz invariant quantities: 
the total amount of entanglement regarded as the sum of these entanglements is not the same in 
different inertial moving frames. In particular, for any value of the entangling parameters, both 
ss and mm-entanglements are attenuated by Lorentz transformations and their parametric rates 
of change with respect to the entanglements observed in a rest frame have the same monotonic 
behavior. However, for indistinguishable (distinguishable) particles, the change in entanglement for 
the momenta is (is not) the same as the change in entanglement for spins. As a consequence, in 
both cases, no entanglement compensation between spin and momentum degrees of freedom occurs. 

PACS numbers: Special Relativity (03.30.+p); Quantum Mechanics (03.65.-w); Entanglement (03.65.Ud). 



I. INTRODUCTION 

It is known in relativistic thermodynamics that the concept of temperature is observer-dependent because radiation 
that is perfectly black-body in a given inertial reference frame is not thermal when viewed from a different moving 
reference frame [HE]- Therefore, since probability distributions can depend on the frames, many other information 
theoretic quantities such as the Shannon entropy can exhibit such dependence. Quantum information theory usually 
involves only a nonrelativistic quantum mechanics. However, relativistic effects are of great importance also in 
quantum information theory j3]. From a practical point of view, relativistic quantum information is important 
since it may provide a useful theoretical platform for several possible applications such as quantum-enhanced global 
positioning [4 , quantum clock synchronization [3] and quantum teleportation [B]. From a foundational viewpoint, 
relativistic extensions of quantum information theory come from quantum cosmology 3 . Specifically, quantum field 
theory in curved spacctime (black hole physics, in particular) present challenges that everybody who upholds the 
principle that "information is physical" [7 should respond to. A first consequence of relativity on quantum theory 
is the existence of an upper bound, the velocity of light, on the speed of propagation of physical effects. A more 
important consequence of relativity is that there is a hierarchy of dynamical variables [5]: primary variables have 
relativistic transformation laws that depend only on the Lorentz transformation matrix that acts on the spacetime 
coordinates. For example, momentum components are primary variables. On the other hand, secondary variables such 
as spin and polarization have transformation laws that depend not only on the Lorentz transformation matrix, but 
also on the momentum of the particle. As a consequence, the reduced density matrix for secondary variables, which 
may be well defined in any coordinate system, has no transformation law relating its components in different Lorentz 
frames. In relativistic quantum information theory, the notion "spin state of a particle" is meaningless if we don't 
specify its complete state, including the momentum variables [S]. It is possible to formally define spin in any Lorentz 
frame, but there is no relationship between the observable expectation values in different Lorentz frames. Stated 
otherwise, the answers to such questions, asked in different Lorentz frames, are not related by any transformation 
group. Under a Lorentz boost, the spin undergoes a Wigner rotation whose direction and magnitude depend on the 
momentum of the particle. Even if the initial state is a direct product of a function of momentum and a function of 
spin, the transformed state is not a direct product. Spin and momentum appear to be entangled. This implies that 
the spin entropy is not a Lorentz scalar and has no invariant meaning in special relativity [5] . 

Entanglement is one of the key features of quantum mechanics and a deep understanding of its properties and 
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implementations is essential not only to fundamentally advance our understanding of how Nature works, but also to 
design more powerful technologies. Entanglement is a property unique to quantum systems. Two systems (microscopic 
particles or even macroscopic bodies) are said to be quantum entangled if they are described by a joint wave function 
that cannot be written as a product of wave functions of each of the subsystems (or, for mixed states, if a density 
matrix cannot be written as a weighted sum of product density matrices). The subsystems can be said not to have a 
state of their own, even though they may be arbitrarily far apart. The entanglement produces correlations between the 
subsystems that go beyond what is classically possible. From a relativistic point of view, entanglement is an observer- 
dependent concept [10] . In particular since Lorentz boosts entangle the spin and momentum degrees of freedom, 
entanglement may be transferred between them. This may be true for single particles [Sj (where entanglement is 
being considered between degrees of freedom belonging to the same particle) , and for pairs [TTHT4] , where the Lorentz 
boost affects the entanglement between spins. Within such quantum relativistic framework, qubits are realized as 
discrete degrees of freedom of particles: a qubit can be either a spin of a massive particle or a polarization of a photon. 
Once the relativistic transformations of the states of massive particles and photons are given, it can be deduced what 
happens to them when described by observers in relative motion. 

In |12j . it was stated that the increase in spin entanglement comes at the expense of a loss of momentum entan- 
glement, since the entanglement between all degrees of freedom (spin and momentum) is constant under Lorentz 
transformations. The misinterpretation of this statement has lead to several misleading remarks in the literature as 
pointed out in |15j . Actually in [15j it was found that the change in entanglement for the momenta is indeed the 
same as the change in entanglement for spins and no entanglement transfer occurs. However, there, it was considered 
a pair of spin-^ indistinguishable particles with maximally entangled momenta and spins in the rest frame. Hence, 
the possibility of entanglement transfer between different degrees of freedom is jeopardized by this assumption (it 
automatically implies that moving from the rest frame, both ss and mm-entanglements will diminish). This led us 
to consider more general scenarios as explained in the following paragraph. 

In this article, we investigate special relativistic effects on the entanglement between either spins or momenta of 
composite quantum system of two spin-^ massive particles, either indistinguishable or distinguishable, in inertial 
reference frames in relative motion. Specifically, for the case of indistinguishable particles, we consider a balanced 
scenario where the momenta of the pair are well-defined but not maximally entangled in the rest frame while the spins 
of the pair are described by a one-parameter (77) family of non-maximally entangled bipartite states. In particular, 
when r\ = (symmetry in the spin-wave function and antisymmetry in the momentum-wave function) or ?y = 1 
(antisymmetry in the spin- wave function and symmetry in the momentum- wave function), we recover the main result 
appeared in [TS]. For the case of distinguishable particles, we take into consideration an unbalanced scenario where 
the momenta of the pair are well-defined and maximally entangled in the rest frame while the spins of the pair are 
described by a one-parameter (£) family of non-maximally entangled bipartite states. Furthermore, in both cases, 
we show that neither the spin-spin (ss) nor the momentum-momentum (mm) entanglements quantified by means of 
Wootters' concurrence are Lorentz invariant quantities: the total amount of entanglement regarded as the sum of 
these entanglements is not the same in different inertial moving frames. In particular, for any value of the entangling 
parameters, both ss and mm-entanglements are attenuated by Lorentz transformations and their parametric rates of 
change with respect to the entanglements observed in a rest frame have the same monotonic behavior. However, for 
indistinguishable (distinguishable) particles, the change in entanglement for the momenta is (is not) the same as the 
change in entanglement for spins. We conclude that in both cases no entanglement compensation between spin and 
momentum degrees of freedom occurs. 

The layout of the article is as follows. In Section II, we present preliminary material on relativistic Lorentz 
transformations for quantum wave-vectors. In Section III, we define the reduced density matrices of the considered 
composite quantum systems of two spin-| particles (either indistinguishable or distinguishable) in a rest frame and 
analyze the action of relativistic Lorentz boosts on them. In Section IV, we present the relativistic effects on the 
entanglement of a composite system of two indistinguishable particles (two electrons with wave- vector Iv&ee)) with 
a balanced (but not maximal) amount of entanglement in the rest frame between momentum and spin. In Section 
V, we analyze the relativistic effects on the entanglement of a composite system of two distinguishable particles (an 
electron and a muon with wave- vector |\P e/J )) in which we consider an unbalanced scenario. Our concluding remarks 
appear in Section VI. 

II. LORENTZ TRANSFORMATIONS FOR QUANTUM WAVE- VECTORS 

In |16j , Wigner showed that quantum states of relativistic particles are given by unitary irreducible representations 
of the Poincare group, the group of translations and Lorentz transformations (boosts) in the Minkowski space. From 
[16) it is also clear that finite dimensional representations of Lorentz boosts are non-unitary. However, special relativity 
requires that the physics of quantum states of relativistic particles should not depend on the arbitrary inertial reference 
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frame from which the states are observed. Thus, we should expect the states to transform unitarily from one inertial 
reference frame to another. The solution to this paradox resides in the fact that in relativistic quantum mechanics, 
the creation and annihilation operators, as well as the associated mode functions for the quantum field that creates a 
given state, transform under Lorentz transformations by local unitary spin-j representations of the three-dimensional 
rotation group. The key ingredient of such relativistic transformations is the Wigner rotation, a rotation in the 
rest frame of the particle that leaves the rest momentum invariant, which restores unitarity in the transformations 
between relativistic single and multi-particle quantum states. Given the explicit expression of the Wigner rotation, 
the transformation properties of entangled quantum states observed from two inertial reference frames moving with 
constant relative velocity can be described. 

An arbitrary transformation L (A, a) of the Poincare group (or, inhomogeneous Lorentz group), where A denotes a 
Lorentz transformation and a a constant vector, relates the coordinates x^ and x'^ of two inertial reference frames S 
and S' [IT], 

ar" — >V = £ (A, a)x" = f A^ + a M . (1) 

The transformation L (A, a) satisfy the composition law, 

£ (A, a) L (A, a) — L (AA, Aa + a) , (2) 

where the bar in ^ is used just to distinguish one Lorentz transformation from the other. Moreover, L (A, a) induces 
a linear unitary transformation £ (A, a) on quantum mechanical wave- vectors |^), 

|*) —►|*') = tf (A, a) |*>, (3) 

where U (A, a) satisfies the following composition rule, 

U (A, a) U (A, a) = U (AA, Aa + a) . (4) 

For the sake of clarity, let us focus on single particle states and for a more general approach we refer to |17j . Assuming 
to consider only transformations with — in the homogeneous Lorentz group, |*') in g for a single particle state 
reads, 



|*) — ► |*') = U (A) |*) d =l f J ^ E ( W ( A > *)) l*A„, *>) , (5) 

where p and a labels the momentum and the spin (or helicity for massless particles) degrees of freedom, respectively. 
The coefficients T>^) a {W (A, p)) provide a representation of the Wigner rotation W (A, p) belonging to the so-called 
Wigner's little group for angular-j particles (for instance, j = \ for spin-i particles). The Wigner's little group 
element W (A, p) is defined as, 

W(A,p) d ^ L- 1 (A P )AL(p), (6) 

where p M = f (p, p°) with, 

_ mv n E mc 

p= — _= andp = - = — (7) 

\fi T W c >/i-(i) a 

and where (Ap) M = f ^PAi (Ap)°J with fj, = 0, 1, 2, 3. The quantity L (p) is the standard Lorentz transformation such 

that p M = L^k u where W — (0, 0, 0, m) is the four-momentum in the rest frame of the particle being considered. 
Stated otherwise, L (p) takes a massive particle from rest to a 4-momentum p. 

For massive (m > 0) spin-^ particle states it turns out that X^ 1 / 2 ) (W (A, p)) in fc>\ reads [TB] . 



2? (1/2) (W (A, p)) = p(-V2) ( L ( Ap )) £,(1/2) ( A ) V (l/2) ^ _ ^ 

Consider an arbitrary boost characterized by the velocity v — ve with e the normal vector in the boost direction. 
Then, omitting tedious technical details that can be found in [131 118) , it turns out that the explicit expression of 
X^ 1 / 2 ) (W(A, p)) reads, 

I?^ 2 ' (W (A, p)) = cos ^ + i (a ■ h) sin ^, (9) 
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where, 



cos 



cosh | cosh | + sinh ^ sinh | (e • p) 
[5 + 5 cosh a cosh 5 + \ (e- p) sinh a sinh (5] 



and, 



(10) 



n sin 



sinh ^ sinh | (e x p) 
[5 + 5 cosh a cosh d + | (e • p) sinh a sinh S~\ 



The quantities a and 5 in (10) and (111 are such that, 

1 



P 



cosh a = — , and, cosh<5 = 

v/l~52 m 



(11) 



(12) 



dcf 



respectively, where j3 = - with c the speed of light. The non-relativistic limit is recovered for /3 = 0. A more detailed 
description of relativistic transformations of either massive or massless arbitrary multi-particle quantum states can 
be found in HU [Tl|. 



III. GENERAL FRAMEWORK: TWO-PARTICLES COMPOSITE QUANTUM SYSTEMS 

In this Section, we study the relativistic action of Lorentz boosts on the density operators of the composite quantum 
systems being considered in a rest frame S. 

We consider two spin-1 particles A and B, either indistinguishable or distinguishable, with positive mass. We focus 
on two-particle wave-vectors that in a given inertial reference frame at rest S may be written as, 

Ws = \^pa, °a\ pb, <y B ) = \pa, <?a; Pb, ob ) = f \pa, Pb) ® Wa, ^b) , (13) 

where p and a denote the momentum and spin degrees of freedom of each particle. Following [15] , we assume that 
the momentum of particle A is concentrated around pa and that of particle B around p B . Thus, the state for the 
momenta of the composite system A + B is a product state that reads, 

\PA, Pb) = \pa) a ®\pb) b - (14) 

Furthermore, \<JAi g b) represents a state for the spins of the two particles. In a different reference frame S' in relative 
motion with respect to the rest frame S, the composite wave-vector \^) s becomes 

l*>5 ^ l*)s' = E vl^\(W(A,p A ))viVll(W(A,p B ))\Ap A ,a> A ;Ap B ,a> B ) 

<*A' °'b 



= \PA, Pb) A ® 2?i V2) (W (A, PA )) V% l2) (W (A, PB )) \a A , a B ) , (15) 

where A is the Lorentz boost between the two reference frames. The transformed state for the momenta of the two 
particles is given by, 

\PA, Pb) \pa, Pb) A = (\pa) a ® I Pb) b ) A = \pa)\ ® \ Pb) b • (16) 

Moreover, the transformed state for the spins reads, 

Wa, °b) S ^ 2?i 1/2) (W(A, pa))V^ /2) (W (A, PB )) \a A) *b) , (17) 

where W (A, p) is the Wigner rotation for the Lorentz transformation A and p^ = (p, m) with m > 0. In the case 
under investigation, X^ 1 / 2 ) (W (A, p)) is a 2 x 2 unitary (rotation) matrix. 

In the rest of the manuscript, we focus our attention on the quantum Lorentz transformation properties of the 
following two-particle quantum mechanical state vector \ty) that, when viewed from the reference rest frame S, reads 

l*)s = I*) = 4? O&x, Pb,) ® |0) s + \PA 2 , PB 2 ) ® \0') s ] , (18) 
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where the state vectors |0) g and \0') s (with |0) s 7^ |0') s , in general) describe the state of spins of the massive fermions. 

We stress that one of the main consequences of the spin-statistics theorem [19] is that the wave-function (both spin 
and space parts) of a system of identical (indistinguishable) half-integer spin particles changes sign when two particles 
are swapped. This can be achieved in two ways: either the space wave-function is symmetric and the spin wave- 
function is antisymmetric, or the space wave-function antisymmetric and the spin wave-function symmetric. Fermions 
are particles whose wave-function is antisymmetric under exchange. Observe that no measurement can distinguish 
between identical particles so we must always have this exchange symmetry in the wave-function of identical particles. 
For distinguishable particles that exhibit different intrinsic physical properties, no exchange symmetry needs to hold. 

From (18) we imply that just two pairs of momentum values (p Al , Pb x ) and (pa 2 , Pb 2 ) are needed to fully char- 



acterize the momentum degrees of freedom of the composite system under investigation. As pointed out earlier, we 
also consider the working hypothesis that such momenta are concentrated closely enough around their distinct values 
so that this allows us to use orthogonal state vectors for distinguishable concentrations. Within such approximated 
scenario, we can employ a single Wigner rotation for each concentration [TS] . 

In a reference frame S' in relative motion with respect to the rest frame S, the quantum Lorentz transformation of 
the wave- vector \^f) reads 



where the transformed momenta state vector reads, 



\PA 1 ,PB 1 ) ® \l) S > + \PA 2 , PB 2 ) ®|2) 5 , 



and the transformed spin states \\) s , and |2) s , are given by, 



|l>s< 



- ID = f x> (1/2) 



A 



(W (A, PAl ))V% /2) (W (A, PBl )) |0) s = V A ( PAl )V B (p fll ) |0) s , 



and, 



|2) s , = |2) d ^V { l /2) (W(A, P a 2 ))V b 1/2) (W(A, PB2 )) |0') s = V A ( p M )V B ( p B . 2 ) |0% ; 



respectively. Substituting (20 1, (21 1 and (22) into (19), the Lorentz-transformed density operator becomes 

dcf 



P 



U{h)prf (A), 



that is, 



\Ap Al , Ap Bl ) (Ap Al , Ap Bl \® |l) (1| + \Ap Al , Ap Bl ) {Ap M , Ap B2 \® |l) (2| + 
+ \Ap M , Ap B2 ) (Ap Ax , Ap Bl \ ® |2) (1| + \Ap A2 , Ap B2 ) (Ap A2 , Ap B2 \ ® |2) (2| 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



where |1) = \1} S , and |2) = \2) s ,. The reduced spin density operator is obtained from p' by tracing over the momentum 
degrees of freedom and it reads 



dcf „ 

= lr 



momentum 



03') = 



1 



(Ap Al , Ap Bl \Ap Al , Ap Bl ) |1) (1| + (Ap Al , Ap Bl \Ap M , Ap B2 ) |l) (2| + 
+ {Ap M , Ap B2 \Ap Al , Ap Bl ) |2) (1| + {Ap A2 , Ap B2 \Ap M , Ap B2 ) |2) (2| 



ID (II 



12} (2|] 



(25) 



with |1) = 1 1)5, and |2) = |2) s ,. Similarly, the reduced momentum density operator is obtained from p' by tracing 
over the spin degrees of freedom and it reads, 



V dcf ry , 

^momentum — lr spin \P ) — 



1 



(1|1) \Ap Al , Ap Bl ) (Ap Ax , Ap Bl \ + (2|1) \Ap Al , Ap Bl ) (Ap A2 , Ap B2 \ + 
+ (1|2) \Ap A2 , Ap B2 ) (Ap Al , Ap Bl \ + (2|2) \Ap M , Ap B2 ) (Ap A2 , Ap B2 \ 



(26) 



The reduced density operators p' spin and p' momcnt um m ^he moving reference frame S' are the analogues of p sp i n and 
Pmomcntum in the rest reference frame S, respectively. 
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FIG - 1: Qpin M = Cmomontum (>?) vs - V, x = V ■ V = (dash) , <p = ^ (thin solid), <p=1 (thick solid). 



IV. INDISTINGUISHABLE PARTICLES 



In this Section, we study the relativistic effects on the entanglement of a composite system of two spin-^ indistin- 
guishable massive particles (two electrons in the state \^f ee ) in the rest frame S) with a balanced (but not maximal) 
amount of entanglement in the rest frame between momentum and spin. We assume that the wave- vector \^s) i n the 
rest frame S is given by, 



|*«) = l*s) = ~\p Ai ,pb x ) [Vv\4>-) + V^v\<t>+)} +^\pa 2 ,pb 2 ) [Vv\<f>-) 

where \4>±) are the maximally entangled Bell-states, 

dcf m) ± lit) 



(27) 



10=1 



V2 




(28) 



and rj is a real parameter such that < rj < 1. Furthermore, we assume that 

pAx = -PA 2 1 Pb x = -Pb 2 and, p Al = p B2 , (29) 
where \pai) and \pa 2 ) are the eigenvectors of ai A ^ with eigenvalues +1 and —1, respectively, while \pbj_) and \pb 2 ) 

( B) 

are the eigenvectors of a z with eigenvalues —1 and +1, respectively. Thus, wc have 

T 



\PA 1 ) = \0) A ^[ Q ), \PA 2 



\1) A = f°\ \p Bl ) = \1) B ee r) and, \p B2 ) = \0) B = Q 



Observe that from (27) and (30), when r\ — and r\ = 1, | tyg) in (27) reduces to 



!*+> = ^|io)l^ 



^|O1)|0+) and, |*_)^-L|1O)|0. 



V2 



|01)|^-), 



(30) 



(31) 



respectively. Using (27) and (30), it follows that the reduced spin density operator in the rest frame S becomes (for 
further details, see Appendix A), 



(S) 

■ 
r spin 



l r 



'4% ® 4*1 + (1 - 2r/) trW 4 B ) + (1 - 2 V ) ® ^ - a™ <8> ^ 



(32) 
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that is, 



(S) def 
^spin — momentum 



(l*X*|) = 5 J 




1 I 1 1 - 2r] 
1 — 277 1 



(33) 











We quantify the entanglement between the spin degrees of freedom of particle A and B by means of Wootter's 
concurrence [20] . 



C (p£>) = max {°. Ai - A 2 - A 3 - A 4 } 



(34) 



( s) 

Denoting p sp ( n = p S pin, the non-negative real numbers {Ai, A2, A3, A4} denote the square roots of the eigenvalues of 
the non-Hermitian matrix p S pin/5spin where p S pm is the "time-reversed" matrix defined as, 

p spin d ^ f ( CT (^) g, p * spin (4-4) ® CT f )) , (35) 

and p* pin is the complex conjugate of p sp i n - It turns out that the concurrence Cfpg^J = (77, </j) reads, 

cS(t?,v) = |(1-2tj)|. (36) 
Following the same line of reasoning, it can also be shown that the reduced momentum density operator is given by, 



1 



momentum 



4x2 ® 4x2 + (2*7 - 1) 4 A) ® + (2^7 - 1) 4 A) ® * 



that is, 



^momentum = Tr spin (|*) (*|) = 





110 1 2t7 - 1 

2 I 2i] - 1 1 





(37) 



(38) 



(s) 

The concurrence of ^momentum m (38 1 becomes, 

^momentum (V, f) = I (1 ~ %V)\ 



(39) 



From ([36]) and ( 39 1 , we conclude that we are considering a balanced but not maximal initial scenario as far as the 
entanglement in spins and momenta concern. 

In a moving reference frame S' where PAu Pa 2 i Pb x an d Pb 2 are along the x-axis and the quantum Lorentz 
transformations are along the y-direction with velocity v, it turns out that the Lorentz-transformed spin density 
operator in (251 reads, 



(*') = I 

Pspin 4 



4x 2 ® 4x2 + (1 - 277) cos 2v*rW ® af> + (1 - 2 V ) cos 2^ ® erf) - ® <rf) 



that is, 





1 I 



Pspin (»7> ¥>) 2 | (1 - 277)C0S2(^ 







(1 - 2?7)cos2^ 

1 




(40) 



(41) 



(*') 



The concurrence C spin (rj, tp) becomes 



4p2(v, <P) = 1(1-2??) cos2(^| 



Similarly, the Lorentz-transformed momentum density operator in ( 26 ) is given by, 
1 r 



{?') 

/^momentum 



4x2 ® 4x2 + (2?? - 1) cos 2ya^ ® 4 S) + (2*7 - 1) cos 2^ ® erf) - ^) ® 



(42) 



(43) 
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FIG. 2: c£2 (?) = Cifo'lntum (?) vs. p for z = <^ £ [0, 1] : q = (dash), r, = | (thin solid), r, = \ (thick solid). 



that is, 



r momentum (»?) <?) 





1 [ 1 (2t7- l)cos2<^ 

2 | (2rj- l)cos2<^ 1 




(44) 



The concurrence C 



(S') 



momentum 



(77, </j) is given by, 



(SO 



momentum 



(77, y,) = 1(1-27?) cos 2<p\ 



Thus, from (42 1 and (45), we uncover that in the case of identical spin-^ particles, 



C, 



spin 



(77, ip)=C 



(s>) 



momentum 



(45) 



(46) 



We stress that when r\ = (symmetry in the spin-wave function and antisymmetry in the momentum-wave function; 
first Eq. in (31 )) or rj = 1 (antisymmetry in the spin-wave function and symmetry in the momentum-wave function; 
second Eq. in (31)), we recover the main result appeared in |15j . Finally, we emphasize that for indistinguishable 
particles, the change in entanglement for the momenta is the same as the change in entanglement for spins (see Eq. 
( 46 1 , Figure 1 and Figure 2) . 



DISTINGUISHABLE PARTICLES 



In this Section, we present a detailed study of relativistic effects on the entanglement of distinguishable particles 
(an electron and a muon in the state l^efi) in the rest frame S) where we consider an unbalanced scenario where 
the momenta of the pair are well-defined and maximally entangled in the rest frame while the spins of the pair are 
described by a one-parameter (£) family of non-maximally entangled bipartite states. 
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A. Spin density operators in inertial reference frames in relative motion 



We assume that the state vector |0) s = |0') s = |0) in (18) that describes the state of spins of the massive fermions 
is given by, 



\o) s = |o) = \* ± (0) = Vi I in) ± a/I lit) 



(47) 



with < £ < 1 and where in the limiting case of £ = \ we recover standard maximally entangled pure Bell states. As 
a side remark, notice that the matrix representation of the density operator to \tp± (£)) reads 



h± (0 = 1^ (0><^± (01 = 









o ±^(i-o £, o 
V o o o o / 



(48) 



In order to make as clear as possible the effect of quantum Lorentz transformation on quantum state vectors, it is 
convenient to re-express p^ ± (£) in terms of Pauli matrices and study the effect of Wigner rotations used to construct 
the Lorentz transformations rotate Pauli's operators. It turns out that in terms of Pauli matrices, p^, ± (£) in (48) may 
be rewritten as (see Appendix A), 



h± (0 



4x } 2 ® i&l + (i - 20 °i A) ® i B A - a - 20 Ail ® ± 2 v ^(i^o ( Ti' 4) ® 

±2 V^T^O^ ® 4^ - 4 A) ® <4 B) 



(49) 



where ^cr^" 4 "* , cr^, cri' 4 ^ and ^cri S ', &y B \ ci" 8 ^ denote the Pauli matrices for the spin of particle A and B, respec- 
tively. 



As a working hypothesis, we assume that, 

PAi = ~PA 2 and p Bl = -Pb 2 with p Al 



-apB 2 



(50) 



where p Al , PA 2 , Vb 1 and ps 2 are along the i-axis and the quantum Lorentz transformations are along the y-direction 
with velocity v. The constant real coefficient a must be such that (/'electron — Vmuon = f- In the case under 
investigation, 



PAi = -aPB 2 => PA! = apB 2 7 {vA 1 )m e v Al = a-y (v B2 )m IM v B2 

that is (setting the speed of light equal to one, c = 1), 

m e v Al _ am^v B2 



After some algebraic manipulations, we obtain 

am tl v B2 



va ± 



2 ml - ml) v% 2 



m,. 



1 



a 2 m 2 _ m 2 



~-VB 2 - 



In the case being considered, following the line of reasoning presented in [13], the Wigner rotation angle ip reads 

sinh a sinh 5 



tan if 



cosh a + cosh S 



where, 



cosh a 



and, cosh 



g _ P° _ 7 (^particle) "T-C 2 



1 - U 2 

particle 



Therefore, 



tan </?i = tan ip2 



sinh a sinh 5i sinh a sinh <5 2 



sinh S\ 



sinh J? 



cosh a + cosh <5i cosh a + cosh 5 2 cosh a + cosh <5j cosh a + cosh <5 2 



(51) 
(52) 

(53) 

(54) 
(55) 

(56) 
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that is, 



sinh S± = sinh 62 



^electron 



V 1 - Electron V 1 ~ V 



2 

muon 



^electron — ^muon; 



(57) 



where v Al = ^electron and Vb 2 = ^muon- Thus, from (53 1 and (57), internal consistency requires 



V Al = a- 



--VB 2 = VB 2 



B2 



(58) 



which is true provided that, 



rn,, 



0,2 m fi — m e 2 



(59) 



In conclusion, we require 



PAi = -Pa 2 and p Bl = -Pb 2 with p Al 



~PB 2 



(60) 



It finally turns out that the Lorentz transformations are defined in terms of the following T> (p) operators, 



W If (IS tp If (r>% 

T> A (p Al ) =cos- -isin-cr^ >, V B (p Bl ) = cos - + 1 sin -a\ >, 



T>A (pa 2 ) = cos I + isin -0"^, X> B (p B2 



¥ Vim 

cos ism —<j\> 

2 2 z ' 



(61) 



where </? is the angle that characterizes the Wigner rotation. From (61), it turns out that the Wigner rotations 
W (p Al ) and W {p B2 ) are rotations by ip around the £-axis and W (paT) and W (p Bl ) are Wigner rotations by — (p 
around the z-axis. Thus, substituting (f6ll) into (21) and (22) and using (47), the Lorentz transformed spin density 



operator p' spin = /5^ ± (£, <p) in (25) becomes (for further details, see Appendix B), 



4± (€, v) = 



4x2 



or, after some algebra, 



4x2 + (1 - 20 ® 4x2 ~ (1 - 20 /£> 2 ® <ri B) ± 2^UT^0 cos 2^ ® <r< B) 
±2^(1 -0 cos 2^4 A) ® a ( v B) - <j (A) ® a (B) . 

f>4,± (0 <fi) = h± (0 cos 2 <p + /5^ T (0 sin 2 <p, 



(62) 
(63) 



with pTj,± (0 given in (49). 



B. Momentum density operators in inertial reference frames in relative motion 

In agreement with |15j . we regard the momentum states \p A , p B ) for the composite quantum system under investi- 
gation as two-qubits momentum states. We use = (ai A \ cr y A \ crSf^ ^ and = (ai B \ &y B \ tfz^ ) to denote 
the vectors of Pauli matrices used to describe the momentum qubit for particles A and B, respectively. Specifically, 
we assume that \p Al ) and \pa 2 ) are the eigenvectors of d z A ^ with eigenvalues +1 and —1, respectively. Thus, we get 



( ,') and. \Pa,) = \1) A =(° 1 ). «i4> 



— * ( B) 

Similarly, we assume that \p Bl ) and \Pb 2 ) are the eigenvectors of a z with eigenvalues +1 and —1, respectively. Thus, 
we obtain 

\P B1 ) = \0) B = Qand, \p B2 ) =\l) B =f°\ (65) 
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Combining Eqs. (18), (47 1, (64 1 and (65), the quantum mechanical wave- vector in the rest frame S becomes, 

1 

The Lorentz transformed quantum states corresponding to \pA t ), \pa 2 ), \pb-l) an d \pb 2 ) are given by, 



l*eM> = \*S] 



1 \PA 1 ,PB 1 )\V^ r l\n)±Vlm) +^=\pa 2 ,pb 2 ) >/T=s\n)±y/i\it) 



S^S' 



\p Al ) A = (cos | -idnfffW) | ) A , |^ 2 ) ^1' |p A2 ) A = (cos | + *sinf |1} 



A - 



(66) 



IPflJ — > IpbJ = (^cos~ +ism-c^ ; J|0) b , |ps 2 ) — > |pb 2 



A / <P . . <P~(B) 

cos — — i sin — <t; y 
2 2 z 



(67) 



Upon computation of the inner products (1|2) and (2|1) and using Eqs. (64), (65) and (67), it can be shown that 

Pmomontum in © reads, 



r momentum & <p) 



I {A) 60 I (B) 
J 2x2 J 2x2 



n /-(A) _ -(B) -(A) _ - 

+ cos 2(p I ai (E> CTi — o-j, ' ® a 
+ (1 - 20 sin 2<p (a£ A) <g> 4 S) + 4 A) ® 



For further details on the derivation of ( 68 ) , we refer to Appendix C 



(68) 



C. Entanglement in inertial reference frames in relative motion 



In what follows, we compute the concurrence of the reduced (mixed) density matrices p' spin (£, <p) in (62) and 

^momentum (£> V) in H' 



i. Concurrence of the Lorentz-transformed spin density matrix 



We quantify the entanglement between the spin degrees of freedom of particle A and B by means of Wootter's 
concurrence [50] • Using (34) and (35), it follows that the concurrence of p^ ± (£) in (48) reads 



(69) 



12 



C(x) 



LO 



09 ■- 



os ■- 



Q.7 




H 1 h 



1 1 1 1- 



0.0 0:1 0.4 0.6 O.B LO 



FIG. 4: <£2 entum K) vs. £, a = e : ip = (dash), <p = ^ (thin solid), ^ = § (thick solid). 



where Cg^ n (0 is the concurrence of the spin density matrix as viewed by an inertial observer located in the reference 
rest frame S. In the moving reference frame S' , it turns out that the matrix representation of p^ ± (£, ip) in (62) reads 



i§± (£, v) = 



o 



\ 

1-f ±V^Tl -t)co$2p 

±^(1-0 cos 2^ £ 

V o o o o / 



and PspinyOgpin with p spin = (£, in ( 70 ) becomes, 



KSpinPspir 



/0 0\ 

£ (1 - Q (l + co s 2 2<p) ±2(1-0 v / C( 1 -0cos2^ 

±2Cv^ (l-0cos2<p £(l-0(l + cos 2 2<p) 

V o o o o / 



The eigenvalues of p ap inp S pm in ( 71 ) are 



Ax ={ (1-0(1 + cos 2<p) 2 , A 2 = £ (1-0(1 - co s 2^) 2 , A 3 = A 4 = 0. 



Finally, the concurrence of the Lorentz-transformed spin density matrix reads, 

c£2 & v) = v / 4Ihw> = \/4e(i - o icos2^i 



(70) 



(71) 

(72) 
(73) 



2. Concurrence of the Lorentz-transformed momentum density matrix 



Following our analysis in the previous Subsection, we quantify the entanglement between the momentum degrees of 
freedom of particle A and B by means of Wootter's concurrence C(/9 momen tum) where ^momentum denotes the reduced 
density matrix obtained by tracing over the spin degrees of freedom. 



In the reference frame S' , it turns out that the matrix representation of /5 momen tum (£; f) m (68) reads 



momentum 



/ 1 cos 2(p - i (1-20 sin 2<p 





\ cos 2ip + i(l - 20 sin 2^9 1 



(74) 
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0.0 0:1 0.4 0.6 O.B LO 



FIG. 5: C^'J (<p) vs. tp for x = <p G [0, 1] : £ = \ (dash), f = | (thin solid), £ = § (thick solid). 



cUld Pmomcntum/^momcntum With Pmomentu 



um /^momentum 



(£, <p) becomes, 



£(£- 1) (1 -cos 2^) + 1 cos2</j-i(l-2£)sin2^ 

l\ 

Pmomentumpmomentum — 2 I Q 

cos2<^ + i(l - 2£)sin2< / j £(£ - 1) (1 - cos2</j) + 1 



(75) 



It turns out that the eigenvalues of PmomentumPmomontum in <|75j) are given by. 
Ai = \\ (l+ Vl-4e(l-Osin 2 2^1 , A 2 



l - (l - \/l-4e(l-C)sin 2 2 r - 



n 2 



, A 3 = A 4 = 0. 



Finally, the concurrence of the Lorentz-transformed momentum density matrix reads, 



a 



momentum 



(£, <p) = Vl-4e(l-C)sin 2 2<p. 



(76) 



(77) 



Let us discuss the main consequences that originate from Eqs. (|73|) and (77). Observe that both the ss and 



mm-entanglements in the moving inertial reference frame S' are functions of either the entangling parameter £ or 
the angle of rotation (p. This implies that the analysis of the monotonic behavior of the entanglements-changes in 
different reference frames is tricky since it involves changes to both £ and ip. From (73 1 it follows that for a given 
value of the rotation angle ip, the ss-entanglement is an increasing function of £ for < £ < | while it is a decreasing 
function of £ for g < £ < 1 (see Figure 3). Furthermore, for fixed values of £, Lorentz boosts with greater angles 
of rotation exhibit an attenuating power on the ss-entanglement that is stronger than that related to boosts with 
smaller angles (see Figure 3). From (77) it follows that, unlike the ss-entanglement, for a given value of the rotation 
angle ip, the mm-entanglement is a decreasing function of £ for < £ < \ while it is an increasing function of £ for 
\ < £ < 1 (see Figure 4). However, like the ss-entanglement, for fixed values of £, Lorentz boosts with greater angles 
of rotation exhibit an attenuating power on the mm-entanglement that is stronger than that related to boosts with 
smaller angles (see Figure 4). For the sake of clarity and without loss of relevant information, we focus our attention 
in the rest of the manuscript on values of the rotation angle ip in the interval [0, 1]. From (73) it also turns out that 
for a given value of the entangling parameter £, the ss-entanglement is a decreasing function of ip for < ip < j while 
it is an increasing function of ip for ^ < ip < 1 (see Figure 5). In particular, unlike the analysis presented in [TS], we 

I I fg'\ 

are able to show from ( 73 ) that the absolute value of the rate of change in ip of Q pin (ip) (that is, its "speed" ) depends 



14 




FIG. 6: cl S 2 cnt u m (<p) vs. <p for x = <p e [0, 1] : £ = § (dash), £ = 1 (thin solid), £ = § (thick solid). 



on the degree of entanglement of the initial state in the rest frame: the spin-entanglement degradation occurs at a 
faster rate for states with a higher degree of entanglement, 



l(& 



spin 



dtp 



> 



d c£2 (y ) 

dp 



£=£r, 



(78) 



with > £ m . Furthermore, for fixed values of <p, to lower values of £ correspond lower values of the ss-entanglement 
(see Figure 5). From (77) it follows that, like for the ss-entanglement, for a given value of the entangling parameter £, 



the mm-entanglement is a decreasing function of <p for < tp < j while it is an increasing function of <p for f < <p < 1 
(see Figure 6). However, unlike the ss-entanglement, for fixed values of ip, to lower values of £ correspond higher 
values of the mm-entanglement (see Figure 6). In particular, unlike the study appeared in [15], we can state from 



(SO 



that the absolute value of the rate of change in ip of C momcntU m (tp) depends on the degree of entanglement of 
the initial state in the rest frame: the momentum-entanglement degradation occurs at a slower rate for states with 
higher degree of entanglement, 



dC 



mo mentum 



< 



dC 



momentum 

d<p 



(79) 



with £m > £m- We also stress that our analysis shows that no sum of ss and mm-entanglements is conserved in any 
pair of inertial reference frames in relative motion for any value of the entangling parameter, 



momentum 



(80) 



Finally, it turns out that both ^-changes and (^-changes of variations of either ss or mm-entanglements, AC sp i n and 
ACmomentum, respectively, exhibit the same monotonic behavior on the permitted range of values for each parameter. 
More explicitly, it follows that 



and, 



sign 



sign 



dAC sl 



dp 



spin 



= sign 



— sign 



9 AC, 



momentum 



dt; 



a ac, 



momentum 



dip 



(81) 



(82) 



15 




FIG. 7: 8AC W*' ^ (dash) and 8AC -°"- K ' ^ (solid) with x = £ vs. $ for (p = f y 



where the variations AC sp i n and A C mom entum 



are defined as, 



and, 



AC, 



momentum 



momentum ^momentum 



^0 



(83) 



(84) 



respectively. For an illustrative justification of Eqs. ([81]) and (|82j), see Figures 7 and 8, respectively. This last finding 
implies that no entanglement compensation between spins and momenta occurs and neither it is required by the 
Lorentz invariance of the joint entanglement of the entire wave-function since no analytical constraint exists between 
this quantity and the ss and mm-entanglements. 



VI. FINAL REMARKS 



In this article, we have analyzed special relativistic effects on the entanglement between either spins or momenta 
of composite quantum systems of two spin-^ massive particles, either distinguishable or indistinguishable, in inertial 
reference frames in relative motion. 

Indistinguishable particles with balanced (but not maximal) entanglement configuration. For the case of indistin- 
guishable particles, we consider a single balanced scenario where the momenta of the pair are well-defined but not 
maximally entangled in the rest frame while the spins of the pair are described by a one-parameter (rf) family of 
entangled bipartite states (see Eq. (27)). We find out that in any scenario neither the ss nor the mm entanglements 
are Lorentz invariant quantities (see Eqs. (42) and (45 1). In particular, for any value of the entangling parameters, 
both ss and mm-entanglements are attenuated by Lorentz transformations. Moreover, the change in entanglement 
for the momenta is the same as the change in entanglement for spins (see Eq. (]46[ ) and Figures 1 and 2). In particular 
(see Eq. (31 )), when r] = or r] = 1, we recover the main result appeared in [15] . 

Distinguishable particles with unbalanced entanglement configuration. For the case of distinguishable particles, 
we consider an unbalanced scenario where the momenta of the pair are well-defined and maximally entangled in 
the rest frame while the spins of the pair are described by a one-parameter (£) family of non-maximally entangled 
bipartite states (see (66 1). We present an extensive investigation of the behavior of two-parameters (the angle of 
rotation tp and the entangling parameter £) ss and mm-entanglements for a two-particle quantum composite system 
of massive fermions in inertial reference frames in relative motion. We show that in any case neither the ss nor the 
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0.0 0.2 0.4 0.S 0. 



eAc Bpin (|, g) 



I \I(t. n: "'"Y^ 1 y; (solid) and aACm ° m °" tum (^- y ) (dag^) w ;th x = <p vs. </? for f = |. 



mm entanglements quantified by means of Wootters' concurrence were Lorentz invariant quantities (see Eqs. (73 1, 



(77)): the total amount of entanglement regarded as the sum of these entanglements is not the same in all inertial 
frames. Furthermore, for any value of the entangling parameter, both ss and mm-entanglements are attenuated by 
Lorentz transformations (see Figures 3, 4, 5, 6) and their parametric rates of change with respect to the entanglements 
observed in a rest frame have the same monotonic behavior (see Figure s 7 a nd 8). In particular, in Appendix D, we 
also consider an additional scenario of distinguishable particles (see Eq. (D3 )) where in the limiting case of £ = \ and 
m e = to m = to (indistinguishable particles, two identical electrons), the main result of |15j is reproduced. However, 



unlike the main finding appeared in |15] and our generalized result in (46) which holds for indistinguishable particles, 
we show that in general the change in entanglement for the momenta is not the same as the change in entanglement 
for the spins when considering distinguishable particles. Furthermore, unlike the main result presented in |12j . we 
provide clear evidences that the rate of entanglement changes in both parameters £ and tp are in the same direction 
(not in the opposite direction). Surprisingly, even allowing for non-maximal entanglement in the spin degrees in the 
rest frame S (while keeping the momentum entanglement maximal in S), no entanglement compensation from the 
momentum to spin entanglement occurs for any possible pair of inertial reference frames one of which at rest. 

In conclusion, although a thorough characterization of relativistic properties of quantum entanglement is far from 
being achieved, our study provides general enough results towards such direction. 
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Appendix A: Decomposition of density operators in terms of Pauli matrices 



In this Appendix, we derive Eq. ( 49 1 . Consider a two-qubit system with Hilbert space H = 



§ C 2 and compu- 
tational basis B-u = {|00) , 1 01} , |10yj~|ll)}. It can be shown that a general two-qubit state p can always be written 
up to local unitary equivalence to a state of the following form [ST] , 



1 



P = 



hx4 + a ■ a (g) I 2 x2 + hx2 <8> b ■ a + <ko~% ® c?i 



(Al) 
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that is, after some algebra, 



/ a>3 + &3 + C3 + 1 61 — ib 2 



bi + ib 2 
a± + ia 2 
V Cx - c 2 



&i — ia 2 Ci — C 2 

a 3 - fo 3 - c 3 4- 1 ci + c 2 ai — «a 2 

Ci + c 2 fe 3 - a 3 - c 3 + 1 61 - «6 2 

ax + ia 2 bi + ib 2 c 3 - b 3 - a 3 + 1 



(A2) 



where a and kl 3 are given by a = ( ai, a 2 , a 3 ), 6 = (61, 6 2 , 6 3 ) and a — (a x , a y , a z ) is the operator vector of Pauli 
matrices. Equating p^, (£) in (48) to (A2|, we obtain 



a 1 =a 2 = b 1 =b 2 = 0, c x = c 2 = 2^(1-0. «3 = 1 - 2£, 6 3 = - (1 - 2f) , c 3 = -1. 
Therefore, in terms of Pauli matrices, p^ + (0 becomes 



1 



4x } 2 ® 4x } 2 + (i - 20 ^ ® 4 2) 2 - (1 - 20 /&> a ® ai 2) + 2yea 3 e)^ 1) ® ai 2) + 



Observe that for £ = I we obtain, 



» 2 



1\ 1 

4 



+2V£(l-£)a 



Similarly, equating p^_ (0 in (48 1 to (A2), we get 

01 = a 2 = h = b 2 = 0, Cl = c 2 = -2^(1-0, a 3 = 1 - 2€, 63 = - (1 - 2£) , c 3 = -1. 
In terms of Pauli matrices, p^_ (0 reads, 



h- (0 = 



r /■(!) ^ r( 2 ) 

J 2x2 



i$ 2 + (1 - 20 4 X) ® 4 2 x 3 2 - (1 - 20 j&> 2 ® ai 2 ' - 2yor^oai 1) ® ai 2 ^ 



-2^ (1-0 



1\ 1 



and for £ = | it becomes, 
In summary, we have shown that 



® 4x } 2 - 4 1} ® 4 2) - ® ^ 2) - ^ ® 4 2) 



h± (0 



1 



rC 1 ) ^ 7-( 2 ) 



4 2) 2 4- (1 - 20 ^ ® 4 2) 2 - (1 - 20 J&> 2 ® ai 2) ± 2y£(1^0 ( Ti 1) ® a< 2) 4- 



±2^(l-04 ) ® * 



(1) ^ „(2) _ fl) - (2) 



Oz ® cri 



and in the limiting case of £ = i we recover the case of maximally entangled pure Bell states. 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



Appendix B: Lorentz transformation of the spin density operator 



In this Appendix, we derive Eq. (62 1. Substituting (21) and (22) into (25), we get 



(BI) 
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where /5 sp i n equals p^± (£) in (49). Using (61), we observe that the following transformation rules hold 



V A ( PAl ) a^v\ ( PAl ) = (cos | - i sin \o^) a™ (cos f + i sin = cos <pa?> + ^ ^ , 



V B (p Bl ) oWvl ( PBl ) = (cos | + zsin ( CO s | - isin f <rf >) = cos^f> - sin^f >, 



( P a 1 )4 A)2? a(pa 1 )= (■ 



cos - / sin — oi A) ) cr, ( , j4) (cos - + isin -a^ 
2 2 / y V 2 2 z 



2 



sin ^(T^ + cos <p<Ty A ^ , 



V B ( PBl )a^V B ( PBl )= (cosl+zsinlaf))^) (cos|-^sin| ( 7( s )) = sin^af) + cos ), (B2) 



and, 



Va {p M )o x A) v\ (pa 2 ) = (cos | + * sin (cos | - isin |a 



cos ipa x A ^ — sin ycr^ A ^ , 



V B (Pb 2 )4 B) V b (pb 2 ) = (cos 



cos <p<r x B ^ + sin ^Cy 5 ^ j 



V A {p M )a^V\{p M )=(cos^ + 



i sin ^a {A) \ a[ A) (cos ^ - z sin ^a^ A A = sin ipa x A) + cos ipa^ A) , 



^b (p Ba ) (Pb 2 ) = (cos I - i sin |af >) ^ (cos 



^ + i sin ^crl B) ) = - sin (pa x B) + cos ^ y B) . 



(B3) 



Therefore, using (B2) and (B3), after some tedious algebra it follows that Pg pin in (Bl) becomes 



spin 



4x2 ® 4x2 + (1 - 20 vi A) ® 4x2 - (1 - 20 4x2 ® ± 2^(1-0 cos2^i A) ® 4 S } + 
±2 v /e(l-Ocos2^4 A) ® 4 B) ~ ^ ® °i B) - 



(B4) 



that is, ^ pin equals pfl (0 in |62j) 



Appendix C: Lorentz transformation of the momentum density operator 



In this Appendix, we derive Eq. (68). Using Eqs. (64), (65) and (67), it follows that 



\Kp Al ) (Ap Al | 



f (A) - (A) 
'2x2 i "2 



|Ap Bl ) (Ap Bl \ 



j(B) -(B) 
J 2x2 + Oz 



|Ap A2 ) (ApA 2 



1 2x2 °z 



-, |ApB 2 )(Api 



J 2x2 °2 



(CI) 



and, 



~(A) . .-(A) -(A) .-(A) 
|A Pj4l > (ApA 2 | = r— ^ , |A P A 2 ) (Ap Al | = 7T^—- 



~(B) . .-(B) -(B) .~(B) 

I a - \ / a - I <?x +UJy | . _ . , . _ | CTi - ICTy 

|Ap Bl > (Ap S2 | = ' ' ^ = 2~ " 



(C2) 
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Using dcTJ and |C2|, /3m 0mentum in d68b becomes 



1 



/^momentum 



AA) . -(A) 

!1|1> I 2X Y^ 



(1|2) 



(2|1) 



(2|2) 



j(^) _ 5 (A) 
'2x2 °z 



J 2X2 °z 



(C3) 



Furthermore, substituting (|47f into (|21J) and p2[ ), it follows that (1|1) = (2|2> = 1 while (1|2) and (2|1) are given by, 
(1|2) = Tr 



/0 0\ 

(l-£)(cos2</? + isin2<£) (1 - £) (cos 2^ + ism2ip) 

±y| (1 — £) (cos 2ip — i sin2cp) £ (cos 2<y£> — i sin 2</?) 

V o o oo/ 



and, 



(211) =Tr 



/ 0\ 

(l-£)(cos2y>-ism2</>) ± y/£ (1 - £) (cos 2y> - i sin 2<^) 

±v/£~ (1 — £) (cos2(/? + isin2<p) £ (cos 2<,£> + i sin 2<^) 

V o o oo/ 



respectively. Thus, substituting (C4) and (C5) into (C3|, we obtain 



cos 2ip — i (2£ — 1) sin 2cp, 

(C4) 

cos 2iy9 + i (2£ — 1) sin 2ip, 

(C5) 



momentum 



4x2+^ 



+ [cos - i (2£ - 1) sin 2ip] 



+ [cos 2<^j + i (2£ - 1) sin 2ip] 



jOt) _-(A) 

'2x2 °z 



'2x2 °2 



(C6) 



After some algebra, /3' momentum may be finally rewritten as p' momelLtum (£, V?) in (68) 



/^momentum — ^momentum 



-^2x2 ® ^2x2 + COS 



„ /-(A) _ -(. 
2<p I <Tx ® &x 

(1 - 2£) sin 2</3 ® + cr 



(C7) 



As a final remark, we stress that the density matrices (D9) and (C7) cannot be parametrized by means of the 
decom 
space 

. 1 



decomposition in (All. For them it is needed the most general parametrization for a two-qubit system with Hilbert 
I 2 (g> C 2 . It can be shown that any state p for such a system may be parametrized as [22], 



9 : 



^4x4 + C ■ a ® I 2 x2 + hx2 ® d ■ 3 + 2J IjkGj ® Cfc 

J", fc=l 



(C8) 



where c, d G K 3 and 7.^ are rectZ numbers. 



Appendix D: Additional example for distinguishable particles 



In this Appendix, we consider the case of distinguishable particles with an unbalanced scenario where the momenta 
of the pair are well-defined and maximally entangled in the rest frame while the spins of the pair are described by a 
one-parameter (£) family of non-maximally entangled bipartite states. Unlike the distinguishable case considered in 
Section IV, we consider a composite quantum system described in the rest frame S by a total wave-vector that differs 
from ( 66 ) and characterized by a set of momentum constraints that differs from that in ( 50 ) . 



Consider the following set of working hypotheses on momentum states, 



pA t = ~Pa 2 , Pb 1 = -Pb 2 and, p Al = ap B2 , 



(Dl) 
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where a = This expression for a is obtained by imposing that the two angles of Wigner's rotations if a and <fs be 
equal (for an explicit derivation of this specific expression for a, we refer to Eq. (59) in Section V). The momentum 
states and \pa 2 ) are the eigenvectors of ai^ with eigenvalues +1 and — 1, respectively, while \pb-i) and \pb 2 ) 

( B) 

are the eigenvectors of a z with eigenvalues —1 and +1, respectively. Thus, we consider 

WaJ = \0) A = Q, \pa 2 ) = \1) A = \p Bl ) = \1) B = (T\ and, \p B2 ) = |0) B = Q. (D2) 

In such case, the wave-vector | v E')c in the rest frame S reads, 



= Ws = ^ [WauPbJ <8 |0) s + |pa 2 , Ps 2 ) ® |0) s ] 



i io) [yr^itDi v^iit) 



i 



Note that for £ = |, in (D3) becomes the wave- vector studied in 



l* ± )H f -L|io}|^ ± ) T -L|oi)|0 ± ) 



(D3) 



(D4) 



with \4>±) defined in (28). Focusing on |\&+) (the same results are obtained using and following the same line of 
reasoning presented in Section V, we obtain that the spin and momentum reduced density operators in the rest frame 
S become, 



p {s) (0 

"spin vs; 



4x2 ® 4x2 + (1 - 20 vi A) ® 4*2 (1 - 20 4x2 ® + 2y£(W)ai A) ® *™ + 
+2^(1-04^®^ 



.(B) 



/ \ 

o i-£ (i - o 

o V£(i-o £ o 

V o o o o / 



(D5) 



and, 



{S) 

r momentum 



1 r 

4 



(A) 



2x2 ^2x2 


















1 


-1 








-1 


1 


















(D6) 



respectively. The concurrence of p sp[n 



(rest-frame) / ^\ i (rcst-framc) 

(?] ana p momcn t um are, 



Cgn (0 = and C, 



(S) 



momentum 



= i, 



^spin 



(0 = 



+2^ (1 - £) cos 2<p<T ( y A) ® 4 W) - a 



r (B) 



(A) _ (B) 



(D7) 



respectively. The Lorentz-transformed spin density operator in (25) reads, 

4x2 ® 4x2 + (1 - 20 ® - (1 - 20 ® > + 2 v /£lT^£) cos 2^ } ® + 



/0 0_ 0\ 

o i -£ a/£(i-0 cos2 <p o 

VC(l-0 co s2^ £ 

V o o o o J 



(D8) 
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Similarly, the Lorentz-transformed momentum density operator in ( 26 1 is given by 

r(A) 



r momentum \S/ ^ 



Ik 



( B) 

^2x2 — cos tyVx 1 ' ® & 



~J A ) to 



COS 2(^(7 



+ (1 - 2£) sin 2tp (a 



~J A ) to _ to 





1 

[cos 2ip — i (1 — 2£) sin 2y>] 





[cos2^ + i(l-2£)sin2^] 

1 




(D9) 



Furthermore, omitting technical details, it can be shown that the concurrences of the Lorentz-transformed spin and 
momentum density matrices in the moving reference frame S' are given by, 



C$}(£, <p) = V / 4ClW)|cos2^| 



and, 



respectively. 



( S >) 



-'momentum 



(£, p) = Vl-4£(1-£W 2< ^' 



We remark that although Eqs. (D3) and (Dl| differ from Eqs. (66) and 



pLref 
nd (|D11 



(D10) 



(Dll) 



respectively, we have obtained the 



same final conclusions in both distinguishable cases considered: Eqs. (D10) and ( |D11 | are equal to Eqs. (73) and (77), 
respectively. We also point out that when £ = I and m e = (two electrons), we recover the main result appeared 
in |15j . Finally, we emphasize that for distinguishable particles, the change in entanglement for the momenta is not 
the same as the change in entanglement for spins (see Eqs. (D10) and (Dll) in addition to Figures 3, 4, 5, 6). 
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